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Abstract 

The KARMEN experiment at the spallation neutron source ISIS used from /x+-decay at rest in 
the search for neutrino oscillations u^—* in the appearance mode, with p ( , e+ ) n as detection 
reaction of Vf,- In total, 15 candidates fulfill all conditions for the v^. signature, in agreement with 
the background expectation of 15.8±0.5 events, yielding no indication for oscillations. A single 
event based likelihood analysis leads to upper limits on the oscillation parameters: sin^(26) < 1.7- 
10-3 for Am2>100eV2and Am2<0.055eV2 for sin2(2e)=l at 90% confidence. Thus, KARMEN 
does not confirm the LSND experiment and restricts significantly its favored parameter region for 

PACS numbers: 14.60.St, 14.60.Pq, 25.30.Pt 
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I. INTRODUCTION 



The study of neutrino masses and mixing originating from extensions of the Standard 
Model (SM) is one of the most interesting issues in particle physics which has also consid- 
erable impact on astrophysical and cosmological problems. For example, neutrino masses 
in the range of a few eV would mean a significant contribution to the matter content in 
the universe. In addition, understanding the mass and mixing scheme of neutrinos is a very 
promising tool to improve our knowledge on mass generating mechanisms for all elementary 
particles. 

A very sensitive way of probing neutrino masses and the mixing between different neutrino 
flavors is the search for neutrino oscillations. The experimental progress in this fleld during 
the recent years has been remarkable, yielding strong evidence for neutrino oscillations 
from investigations of solar and atmospheric neutrinos. The long-standing problem of the 
solar i/-deflcit, observed by different experiments [|l| including the latest results from the 
Sudbury Neutrino Observatory (SNO) 0, is consistently explained as the transition of z/g 
into other active neutrino flavors 0,0. In addition, the atmospheric neutrino anomaly gives 
evidence for neutrino oscillations, namely for I'^^Vx disappearance oscillations p. Due to 
the precision measurements of the Super-Kamiokande experiment, the oscillation channel 
z/^— > Vj. is strongly favored Q]. 

Despite the convincing results from solar and atmospheric z/-oscillation experiments, all 
indications for oscillations are obtained by searches in the disappearance mode. Up to now, 
there is only one piece of evidence for z/-oscillations in the appearance mode: the LSND 
(Liquid Scintillator Neutrino Detector) experiment at the Los Alamos Neutron Science 
Center (LANSCE) reported 1995 initial results of the search for z?^— > Ve oscillations with 
produced in /x+ decays at rest §]. Supported by a positive signal in the channel 0, 

updates with increased statistics [[l^, |Tl[| underlined the evidence of an observed excess 
but also reduced the original signal strength. The Ve signal is explained as originating from 
z/^^ z/e oscillations with an oscillation probability P = (0.264 ± 0.067 ± 0.045)% [0]. 

Due to the sensitivity region of LSND, these flndings suggest rather high mass differences 
of Am^ > 0.1 eV^, which would imply significant contributions of neutrinos to the cosmolog- 
ical problem of dark matter. Due to the high Am^ scale it is not possible to accommodate 
all three evidences (solar, atmospheric, LSND) with their distinct regions of Am^ within 
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the framework of the SM with its three neutrino flavors, extended by allowing for non-zero 
neutrino masses. Proposed solutions to this problem include e.g. the incorporation of a 
sterile neutrino state []13|, |1^, 0, supersymmetry |T6| or CPT violation |T^. These deep 



impacts on particle and astrophysical aspects therefore require a thorough and independent 
test of the t'^^ Ve evidence of LSND. 

This paper describes the search for v^—^ Ve oscillations by the KARMEN (KArlsruhe 
Rutherford Medium Energy Neutrino) experiment, which was located at the highly pulsed 
spallation neutron source ISIS of the Rutherford Laboratory (UK). The results presented 
here are based on the final data set recorded with the full experimental setup of KARMEN 2 
from February 1997 until March 2001. 

The KARMEN experiment took data, in a different experimental configuration (KAR- 
MEN 1), since 1990. In this first period, the data analysis was focused on the investigation of 
neutrino-nucleus interactions |jl8|, |19|, [2^] , but also on the search for the oscillation channels 
y^^Vf, [0] and Vf,^Vx Other searches of non standard model physics such as new 



particles in pion decay [Q, lepton flavor violating pion and muon decays p4[ or non V-A 
contributions to the muon decay — ^ e^ + z/g + z/^j |2^ were also performed. Here, we 



report on the most sensitive channel, the search for v^-^ z/g oscillations. 

The paper is organized as follows: Section II describes the neutrino source ISIS and the 
KARMEN detector, after which, in section pTl] , the v^-^ Ve oscillation signature is presented. 
Section ^ defines some general event requirements for the identification of z/g-induced events 
in the data analysis. We discuss the background in section 0. The final event sample 
together with the final data cuts and background expectations is given in section The 
data analysis is described in detail in section [VII| together with the presentation of the final 



^[i-^ results. A detailed discussion of the results with respect to the LSND evidence and 
the negative results from other experiments follows in section [VII C . 



II. EXPERIMENTAL CONFIGURATION 
A. The neutrino source ISIS 

The pulsed spallation neutron source ISIS of the Rutherford Appleton Laboratory uses a 
rapid cycle synchrotron to accelerate protons up to 800 MeV with a design beam current of 
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FIG. 1: (a) Time and (b) energy distribution of neutrinos at the ISIS beam stop for a beam 
current of / = 200 /iA: z^^ from /i^ decay (solid), i>e from ^~ decay (dashed). 



/ = 200 /iA. The protons are extracted from the synchrotron with a repetition frequency of 
50 Hz as a double pulse, consisting of two parabolic pulses, with a width of 100 ns and being 
separated by 325 ns in time. When the 800 MeV protons hit the water cooled Ta-D20 target 
(0.0448 ± 0.0030) 7r+ per incident proton are produced p6|. Production of three distinct 
neutrino flavours z/^, and occurs via the tt^- //^ decay chain in the beam stop: 



26 ns 



2.2 fis 



The 7r+ and fi^ are stopped within the heavy target and decay at rest. The unique time 
structure of the ISIS proton pulse allows a clear separation of z/^ induced events from z/^ 
and z/e induced events. Due to the short life time of 7r+ (r = 26 ns) the z/^ production 
closely follows the ISIS proton beam profile. One therefore expects two z/^j bursts within 
the first 600 ns after the extraction of the proton beam. The 2-body decay at rest of 7r+ 
leads to monoenergetic with an energy of Eu = 29.8 MeV. Studies of these are 



published in pO|- the other hand the z/^ and z/e from /i"*" decay are expected to emerge 
on a time scale of a few /xs due to the /i"*" life time of r = 2.2 fis. The time spectrum of 
z/^ and z/g induced events [see Fig. |I] (a)] reflects the Ufe time of /i"*" and thus contains 
additional information to discriminate in the data analysis versus background reactions. 
The z>^ and z^e from muon decay have continuous energy spectra (see Fig.[^). The energy 



spectra are well defined and can be calculated precisely because of the decay at rest 
kinematics and the simple V-A structure of the /i"*" decay. From the three neutrino flavours, 
which are produced with equal intensity and emitted isotropically, the highest mean energy 
is obtained by the z/^, which have the maximum intensity at the endpoint energy of 52.8 MeV. 



The intrinsic contamination of the ISIS z/-beam with v^, is very small. The suppression 
of i^e^pro duct ion follows from the following factors: The stopping of 800 MeV protons in 
the Ta-D20 target produces less vr" than 7r+ (7r^/7r+=0.56). While vr^, which are stopped 
quickly (< Ins), mainly undergo nuclear capture, it is only a fraction of 1.2% which decay 
in flight and therefore become of relevance for the z/g contamination. The following /i~ decay 
at rest in the target station again is suppressed by the efficient muon capture (93% of /i~ 
produced) on the high Z material of the spallation target. This 7r~-/i~ decay chain leads to 
a very small contamination of z/g/ z/^j = 6.4 • 10"^ p6[ with the distributions for z/g in energy 
and time shown as dashed lines in Fig. |l]. The intrinsic z/g contamination is discussed in 



more detail in section |VB3| . The small z/g component in the ISIS z/-beam together with 
the unique time structure of the proton beam allows a high sensitivity search for z/^— > z/g 
oscillations. 



B. The KARMEN detector 



The KARMEN detector is a segmented high resolution liquid scintillation calorime- 
ter, located at a mean distance of 17.7 m from the ISIS target at an angle 100 degrees 
relative to the proton beam. The liquid scintillator is enclosed by a multilayer active veto 
system and a 7000 t steel shielding (see Fig. |^). The hydrocarbon acts as active target 
for neutrino-nucleus reactions (^^C,^'^C,^H). The 65 m^ of liquid scintillator consisted of a 
mixture of paraffin oil (75%vol.), pseudocumene (25%vol.) and 2 g/1 of the scintillating 
additive l-phenyl-3-mesityl-2-pyrazoline (PMP). 

The liquid scintillator volume is optically separated into independent modules by an optical 
segmentation of double Incite sheets. A small air gap between the double Incite sheets of the 
segmentation causes optical total reflection and thus a very efficient transport of scintillation 
light to the ends of the modules, where the scintillation light is read out by a pair of (3" 
VALVO XP 3462) photomultiplier tubes (PMT). Furthermore, gadolinium coated paper has 
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FIG. 2: (a) Front view of the KARMEN detector with details of the central detector region and 
a single module, (b) Side view, the ISIS target is located to the right. 



been put between the acrylic walls for an efficient detection of thermal neutrons. 
Tlie segmentation consisted of 608 modules in total, which are placed inside a rectangu- 
lar tank with the dimensions of 3.53 m x 3.20 m x 5.96 m in length, width and height. 
The central detector consists of the inner 512 modules (each with the dimensions of 
353 cm X 17.7 cm x 18.1 cm in length, width and height), arranged in 32 rows and 16 
columns. A surrounding layer of modules with half the cross section of a central detector 
module defines the inner anti counter. An inner passive shielding of 18 cm thick steel slabs 
surrounds the scintillator tank providing passive shielding and mechanical stability. The 
second layer of active shielding (inner veto) consists of 136 plastic scintillator bars (NEllO) 
with thicknes of 3 cm and lengths ranging from 2.4 m to 3.1 m, which are mounted onto the 
passive shielding on all sides but the bottom side. 

The surrounding steel shielding is built in a modular way out of layers of steel slabs. This 
structure of layers allowed the integration of an outer veto system inside the steel shielding. 
In total, 136 bars of plastic scintillator (Bicron BC412) have been used for the outer veto 
system, which provided also active shielding under the detector. 

This additional outer veto system was installed in 1996, marking the beginning of the KAR- 
MEN 2 experiment. The upgrade of the experimental configuration reduced considerably 
the background level for the Ve search, as it will be outlined in section V A . 

The KARMEN detector is a liquid scintillator calorimeter, optimized for high energy 
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resolution of as = 11.5%/ a/ -E'(MeV). An event information comprises the energy, time 
and position information, as well as the number of addressed modules and their relative 
time differences. A scintillator module hit is accepted if there is a coincidence of signals 
of the photomultipliers at both ends within a coincidence time of ATci = 190 ns (first 
level trigger). The position of the event along the module axis (x-direction) is obtained 
by the time difference between the signals, whereas the energy information is derived from 
the integrated PMT pulses. The absolute energy calibration of the detector is fixed by the 
analysis of the Michel energy spectrum of electrons from the decay of stopped cosmic ray 
muons. The energy calibration is performed for each single module and takes into account 
the individual light output curves of the modules. Module hits within a coincidence time 
AT(7 < 90 ns are combined to one event. Analysis of throughgoing muons allow to calibrate 
the relative times of module hits trei. with an accuracy of 6trei. = 0.8 ns (FWHM). In 
the case of events with more then one module hit, the 3-dimensional position information 
(x,y,z) corresponding to module axis, row and column is constructed by the energy weighted 
average of the single module information. Finally, the event time t relative to the ISIS 
proton beam is recorded. Individual KARMEN modules are synchronized to the ISIS beam 
with an accuracy of 6t < 2 ns, allowing to exploit the ISIS time structure in detail. A beam 
reference time of t=0 is attributed to the time, when the first neutrino enters the KARMEN 



detector. A full description of the detector energy and timing calibration is given in p8 



III. OSCILLATION SIGNATURE 

Neutrino flavor oscillations occur, if the weak interaction eigenstates u^, and Vr are a 
superposition of the non-degenerate mass eigenstates z/i, z/2 and z/3. As the mass eigenstates 
propagate differently, there is a non-zero probability that a neutrino flavor produced via the 
weak interaction (e.g. v^) is detected as another neutrino flavor (e.g. z/g) after a traveling 
distance L. In general, the formalism of the mixing of three flavor and mass eigenstates 
requires a unitary 3x3 mixing matrix f/, often referred to as the Maki-Nakagawa-Sakata 
| 29| matrix Umns- However, the current results in the fleld of neutrino oscillations suggest 



a one-mass-scale dominance dm? = Am^i2 <^ Am^is and Am^ia ^ Am^23 = Am^ with 
Am'^ij = \m\ -m^j\;i,j = 1, ... ,3 fSg, 0, 0, Possible mixing to sterile 



neutrinos as suggested by [O |T1|, W^ is ignored whereas CP conservation is assumed, as 



8 



we shall do in the following. In this case, and since the KARMEN experiment with its 
distance between neutrino source and detection point of L ^ 17 m is a typical short baseline 
oscillation experiment, it is sufficient to simplify the mixing scheme to a 2 x 2 mixing. In 
such a two flavor mixing scheme, the probability P to detect a in an initially pure 
beam with energy E (in MeV) after a path length L (in meters) can be described as: 



P{v,-^ v.)=A. sin^ j (1) 

In a short baseline regime (1/Am^ L/ E <^ l/5m?), contributions to the oscillation 
probability P due to the smaller difference of the squared //-masses, 5m, can be neglected. 
The oscillation amplitude A in (P is a function of the elements of the mixing matrix ?7mns ■ 
For simplicity, we define 



A = sm\2e) (2) 

keeping in mind, that for a comparison of oscillation searches in a different mode than 
i^ti^ appearance, one has to calculate A as the complete function of the 3x3 mixing 
matrix elements. For a review on neutrino masses and mixing and a complete formalism of 
neutrino oscillations see p7|. 



A. Vf. absorption on protons 

Appearance of z/g from Vy,-^ z/g flavor oscillations is detected by the classical inverse beta- 
decay on the free protons of the scintillator: 

z>e + P — n +@ Q=- 1.804MeV 

t 

^therm 

+Gd-^Gd+n@ (n) = 3 

The z/e signature is therefore a spatially correlated delayed coincidence between a prompt 
positron and a delayed 7 event from a (n, 7) neutron capture reaction. 
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1. Positron signal 

For different sets of parameters sin^(2B) and Am^ the oscillation probability P{p^—^ v^.) is 
calculated varying energies and flightpaths. These Vf. energy spectra are then transformed 
into positron spectra by means of the calculated energy dependence of the p ( i^e , e+ ) n cross 



section. The calculation used [39| takes into account weak magnetism and recoil effects, 



yielding a flux averaged cross section of atot = 93.5 x 10^^^ cm^ for the z/^ spectrum from 
/i"'"-decay at rest. Due to the short baseline of (L)=17.7 m, the strongest z/^— >z/e signal 
is expected at Am^ = 2.8 eV^. Figure ||(a) shows the dependence of expected e"*" energy 
spectra for three mass difference values (Am^ = 1, 10, 100 eV^), illustrating the modification 
of the energy spectrum due to oscillation effects. The spectra include experimental response 
functions such as energy and spatial resolutions, threshold efficiencies as well the integration 
of the oscillation probability over the detector volume. The visible energies of positrons 
extend up to 50 MeV with the oscillation signal mostly above 20 MeV. Figure |^(a) also 
demonstrates the power of the detector to discriminate between different values of Am^ in 
case of a positive oscillation signal. 

Apart from the well defined energy spectrum, the time spectrum of e"*" [see Fig. 0(b)], 
resulting from the unique ISIS time structure, discriminates against beam uncorrelated back- 
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FIG. 4: (a) Energy and (b) time distribution of neutron capture events. The energy signal (ex- 
perimental data points) is the sum of p ( n,7 ) d (MC dotted line) and Gd ( n,7 ) (MC dashed line) 
capture. The time between neutron production and capture is quasi-exponential with a time con- 
stant of r ~ 120 /is well reproduced by MC. 

ground. The time distribution of the positrons follows the 2.2 /is exponential decrease of the 
/i^ decay at rest. The positrons are therefore expected in a narrow time window of several 
/is after beam-on-target . 

2. Neutron capture signal 

The delayed event of the induced delayed coincidence arises from one of two different 
neutron capture reactions. Neutrons from p ( z/g , ) n reactions have kinematic energies 
up to 5 MeV and are quickly thermalized. After thermalization, neutrons are captured either 
on protons of the scintillator p ( n,7 ) d or on gadolinium Gd ( n,7 ) , which is contained inside 
the walls of the segmentation. In the first case, a single mono-energetic 2.2 MeV gamma is 
produced, in the latter case, a complex gamma cascade is initiated with a sum energy of 



Neutron capture reactions are monitored in situ during the measurements by investigating 



Y.E^ = 7.9 MeV H [see Fig. |(a)]. 
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the capture reaction 



of stopped cosmic ray muons. This reaction produces neutrons with kinetic energies in the 
few MeV range [^], comparable to the energy of neutrons from the p ( i^e ,e+ ) n process. 
Figure ^(a) shows the measured spectrum of visible energies following a stopped muon in a 
coincidence volume of ~ 1 m^ d'^^;! < 60 cm, |Arow|, |Acol| < 2.5) around the endpoint 
of the muon track. The p ( n,7 ) d peak can be clearly separated from the broad distribution 
of Gd ( n,7 ) signals. The Gd ( n,7 ) signal does not peak at Eq = 7.9 MeV due to the calori- 
metric properties of the single modules. If the 7's from the cascade are spread over different 
modules, missing visible energy can occur due to the thresholds of individual modules. 
The neutron thermalization and capture followed by 7 emission is simulated using the 



GEANT/GCALOR program |42, 43]. The simulated spectra shown in Fig. ^(a) include 
detector response functions and have been adjusted separately to the measured distribution. 
For visible energies below 3-4 MeV the energy resolution, as well as hardware thresholds 
together with the complex topology of a multi-7 event lead to difficulties in describing the 
spectral shape by Monte Carlo simulations. However, since fi~ capture reactions (Eq. ^ 
are measured, the spectral shape of neutron capture events and the total neutron detection 
efficiency can be reliably measured, in order to be used for the z/^^ Ue search. 
The experimental as well as the MC generated time difference between the prompt cosmic 
muon and the 7's from the neutron capture is shown in Fig. §(b). The distribution can be 
approximated by a single time constant of r ^ 120 /is, reflecting the thermalization and dif- 
fusion processes of the neutron and the subsequent two competing capture processes. There 
is a slightly enhanced occurence of 7's within the first fis is due to a higher rate of Gd ( n,7 ) 
capture. This is explained by the almost immediate capture of neutrons being produced 
near the walls containing Gd. 



3. Neutron detection efficiency 

The neutron detection efficiency has to be determined accurately in order to calculate 
the expected number of (e+,n) sequences from v^^Ve oscillations. The efficiency is 
determined by monitoring the nuclear capture reactions of stopped muons (Eq. It is 
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given by the ratio of detected neutrons Nn to the total number of produced neutrons 
The number of detected neutrons Nn is given by the number of delayed coincidences occuring 
after a stopped muon. According to the expected neutron capture signal, we require the 
delayed event to occur within a coincidence time 5 < At < 300 /xs with energies -E^ei. ^ 
8 MeV and within a coincidence volume of Vc = 1.3 m^. 

In order to derive the total number of produced neutrons M„, the number N^- of stopped 
the /i~ capture rate Ac, and the neutron emission multiplicity (x) must be known. As 
the charge of stopped cosmic muons cannot be determined for individual tracks, the decay 
time spectrum has been analyzed to derive the charge ratio /i^/ = Rf^ = 1.28 ± 0.03 and 
thus the number A*"^- of stopped fi~ is known from the measured number of stopped muons 
Nf,. With a total fi' capture rate of A*°* = (38.4 ± 0.4) • lO^^s"^ on g corrected for 
the abundance of ^^C and ^^O in the scintillator, an average probability per stopped fi~ of 
a" = (64.1 ± 1.3) ■ 10^^ is derived for processes with neutron production. 
The derived neutron detection efficiency e from these values 

Nn-{1 + R,) 

must then be modified in two aspects: 

(1) Due to multiple neutron emission (x) = 1.07 (see Eq. ^), the derived efficiency emust 
be corrected to the single neutron expectation from the p ( i^e , ) n reaction. 

(2) As the identification of the muon stop point can lead to ambiguities for tracks, which 
stop close to the borders of the detector, a restricted fiducial volume of the detector to 
the stop points of muons (Ix^topl < 150cm, the outermost module layer removed) is ap- 
plied. The detection efficiency e is then extrapolated to the entire detector volume using 
GEANT/GCALOR simulations. 

A complete description of the analysis of muon capture reactions with the KARMEN 



detector and the derivation of the neutron detection efficiency is given in [45 



Taking all effects into account, the neutron detection efficiency En amounts to: 

En = 0.42 ± 0.03 (5) 

This value is the neutrino flux weighted average of the entire KARMEN 2 measuring period 
as shown in Fig. |^. 
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FIG. 5: Measured single neutron detection efficiency as a function of time during data taking. The 
liorizontal bars indicate ISIS beam-on intervals, the dotted line shows the neutrino-flux weighted 
average of the neutron detection efficiency, the dashed lines the total systematic error band. 

B. Vf, absorption on carbon 

A second z/g detection reaction is the inverse beta decay of carbon ^^C ( z/g , e"*" n ) ^^B with 
a Q- value of 16.7 MeV. This detection reaction has a smaller flux- averaged cross section 
^ than p ( z/g , e"*" ) n . In addition, the number of target atoms Nt in the scintillator is 



smaller than the number of free protons (see table 1). It thus contributes about 5% to the 
detection of Ve- The GEANT 3.21 Monte Carlo simulation of ^^C ( z/g , e+ n) ^^B is included 
in the total number and spectral shape of expected (e+,n) sequences from v^-^ Ve oscillations 
[Fig. 3(a)]. 



p ( z^g , e+ ) n 
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C ( Pg , e+ n) ^^B 



(j{i>e contamination) 



4.5 • 10^° 
93.5 • 10-^2 gj^2 
72.0 • 10-^2 



2.5 • 10^° 
8.5 • 10-^2 pj^2 
7.4 • 10-^2 



TABLE I: Comparison of flux averaged cross sections a and target nuclei Nt for detection of fg 
from different sources. 
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IV. GENERAL EVENT REQUIREMENTS 



The special feature of the z/^— Ve signature is its delayed coincidence nature of a prompt 
high energetic positron, followed by a low energetic signal from neutron capture. Before 
enforcing stringent cuts, which correspond to the delayed coincidence nature of the v^. 
detection reaction, we apply loose cuts to the data set, which do not cut into the signal 
region but which strongly suppress background. 

(1) Only sequences of two events are accepted. 

(2) A sequence accepted for further evaluation in the software analysis consists of a prompt 
event and a delayed event which shows the typical characteristics of neutron capture 
events. In particular, this means that the delayed signal occurs within At < 500 yus after 
the prompt event and has energies less than Eii^i, < 8 MeV. A coincidence volume of 

= 1.3 m'^ is required. 

(3) Neither the prompt event nor the delayed event must have any hits in the multilayer 
veto system. 

(4) The prompt event must have energy Ep^, > 11 MeV. 

(5) There must be no activity in the detector system preceding a prompt event. The 
history of all activities in the detector system (total trigger rate Ttot ~ 13 kHz) are stored 
by a time stamp and a bit pattern word, which allows the decryption of addressed detector 
parts. Requesting no activities preceding an event in the main detector, inner veto or inner 
anti counter in the previous 24 /is (14 /is for the outer veto system) eliminates most of the 
cosmic induced background with short time correlations, as shown in figure |^. 

(6) There must be no stopped cosmic ray muons in the central detector preceding a 
prompt event. With a rate of ~ 160 Hz the hardware trigger identifies stopped muon 
in the central detector. A 10 /is hardware dead time is then applied and the event time 
and stopping position of the muon are stored, thus providing information for the offline 
analysis to detect spatial correlations between an event and preceding stopped muons. 
Prompt events of a potential Ue coincidence are rejected, if they occur within At < 40 /is 
after stopped muons anywhere in the central detector, after up to At < 500 /is within 
a coincidence volume of Vc = 1.3 m'^ (/i~ capture with n-emission), or if they occur in a 
coincidence volume of Vc- = 0.5 m'^ for time differences At < 100 ms (yU~ capture with 
subsequent ""^^B /3 decay). 
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FIG. 6: Rate of events following in a time difference At to the last preceding event (a) in main 
detector, inner veto or anti counter (b) in outer veto system. The count rate suppression for time 
differences At < 15 /xs in Fig. (a) is caused by hardware and software deadtimes as well as read-out 
dead times. 



(7) In the case of events with more than one addressed module in the central detector, 
the maximum time difference between the module hits must not exceed ATcmod = 50 ns, 
ensuring that the module hits belong to the same physical event. 
(8) Not more than 10 modules of the central detector must be addressed. 



V. BACKGROUND REACTIONS 

Evidence for flavor oscillations z/^— Ve in the appearance mode requires statistically sig- 
nificant detection of z/g in the time window of in excess of any inherent background. 
While for maximal mixing one expects several thousand oscillation events, a mixing ampli- 
tude 10"'^ < A < 10"^ (as suggested by LSND) could reduce this number to about 10 events. 
Despite the clear oscillation signal and the small ISIS duty cycle, the clear and unambiguous 
detection of such rare Ve events requires a very efficient detection and suppression of the large 
amount of cosmic induced reactions. Benefiting from the three-fold active veto system the 
cosmic background can be suppressed to a level well below the expected oscillation events. 

However, neutrino induced reactions can also induce a background rate. In particular, u^. 
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induced charged and neutral current reactions constitute the largest background reactions 
in the search for z/^— > Ve oscillations. This section discusses both background reactions in 
the z^e search, induced by cosmic rays as well as by neutrinos. 

A. Cosmic induced background 

The cosmic ray induced background reactions are measured in the long beam-off time 
window between the i^-pulses. Taking into account the trigger structure of the experiment, 
which also allows for calibration measurements, the effective statistics for cosmic induced 
reactions in the beam-off time interval is 140 times larger than the narrow time interval for 
the z/-pulse. This factor allows to extrapolate the determined cosmic induced background 
rate with a statistical accuracy of 5% of the neutrino analysis. 

The 7000 t steel shielding of the detector absorbs both the hadronic and electromagnetic 
component of cosmic rays. It is therefore only the muonic component, which can induce 
z/g-like background processes. 

1. Throughgoing muons 

The KARMEN central detector was exposed to a rate of 1.1 kHz of throughgoing muons. 
These muons were detected in the central detector modules, as well as by the active veto 
system. The veto system inefficiency is estimated to be less than 2.2 x 10~^. Delayed 
activities following cosmic ray muons by spallation processes of high energetic muons on 
^^C, are highly suppressed due to the general event requirement 5 (see section |IV|) and can 
be neglected in the z/g search. 

2. Stopped muons 

Stopped muons in the central detector can cause spatially correlated events on the time 
scale of a few microseconds up to several milliseconds. Whereas all /x"*" stopping in the 
detector will decay, a fraction of ac = 7.8% of the stopped fi~ undergo nuclear capture 
reactions in the scintillator. The muon decay produces a spatially correlated electron or 
positron with an energy up to Eq = 52.8 MeV. The time correlation is defined by the 
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lifetime of /i"*" (r = 2.197 fis) and (t = 2.026 /is). With a branching ratio of F - = 0.82, 
the nuclear capture reactions involve neutron production: 

The neutrons are detected by the typical neutron capture events of p(n,7) or Gd(n,7) with 
Eq = 8 MeV and Tcapture ~ 120 /is. This process leads to a contribution to the cosmic 
induced background in the search, which arises from unvetoed muons with short track 
lengths, stopping in the central detector and depositing less than 51 MeV. 
Long lived background arises from muon capture reactions of /i~ 

to the ^^B ground state or 7-unstable levels, through the subsequent /3-decay: 

i^Bg.,. + e- + z7e (8) 

with r = 29.1 ms and an endpoint energy of Eq = 13.3 MeV for the beta-electron. Hence, 
this reaction has only a small overlap in its signature to z/g induced coincidences. Never- 
theless, each event arising in the main detector is checked for preceding stopped muons for 
time differences up to At < 100 ms (general event requirement no. 6) to suppress the beta 
decay, whose electrons otherwise give rise to random coincidences. 



3. Muons near the central detector 



The dominant cosmic ray induced background is due to muon interactions in the 7000 t 
steel shielding blockhouse, which generate highly energetic neutrons. Two different reaction 
mechanisms can be distinguished: 

• fi~ capture on ^^Fe: 

/i- +^^Fe — > ^^-^Mn + x-n + (9) 
Negative charged muons stopped in iron are predominantly captured with a capture 



rate of Ac = (4.411 ± 0.026) x 10^/s The energy transferred to the nucleus in the 



process is between 15 and 20 MeV and therefore above the neutron emission threshold. 
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FIG. 7: Energy distribution of prompt events of cosmic induced sequences. Measurement ignoring 
information (open dots) and using information (full dots) of the outer veto system. See text for 
details on the exponential fits. 

• Deep inelastic scattering (DIS) of muons on ^^Fe: 



Virtual photons radiated from the cosmic muons interact with the iron nuclei and can 
produce spallation neutrons with energies up to a few GeV. On average, 3-4 secondary 
particles with energies above 10 MeV are produced, primarily neutrons and protons. 

Neutrons from deep inelastic scattering can penetrate into the liquid scintillator, causing 
signals with visible energies up to 200 MeV through elastic n-p scattering. After thermal- 
ization the neutrons are captured either on protons or on the gadolinium, yielding capture 7 
spectra, as shown in Fig. ^. Thus, the highly energetic neutrons cause delayed coincidences, 
which are nearly identical to the signature of Ue, as the KARMEN detector has no particle 
identification and cannot distinguish between cosmic induced n-p recoil events and positrons 
from p ( z/g , ) n . The crucial identification of the highly energetic neutrons is achieved 
by the third veto counter system, which is placed inside the steel shielding. Figure |^ shows 
the spectrum of the visible energies of the prompt events, covering the entire energy interval 
of a potential oscillation signal. The delayed events of these sequences follow the expected 
distributions for neutron capture (see fig. ^). 




X + y ■ n + 



(10) 
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FIG. 8: Distribution of time difference 6t between liits in the outer veto and subsequent hits in the 
central detector of cosmic induced background. The Monte Carlo simulation (solid line) consists 
of three components: (i) fast neutrons from DIS, (ii) neutrons from muon capture on iron, (iii) 
stopped muons. 

Figure |^ shows the identification of the processes involved by the time correlation of 
prompt muons and the proton recoil event. The time distribution is measured by the time 
difference St between the hit in the outer veto system caused by the muon and the subse- 
quent hit in the central detector caused by the proton recoil from highly energetic neutron 
interaction. The time distribution shows three components: 

(i) The dominant gaussian shaped distribution peaking at a time difference of 6t = 25 ns 
with an additional enhanced tail distribution, which can be attributed to highly energetic 
neutrons from deep inelastic muon scattering on iron. The time difference for these events 
is equivalent to the time of flight of the neutrons from their point of production in the steel 
shielding to their n-p interaction in the central detector. 

(ii) For time differences 6t > 60 ns neutrons from fi~ stopping in iron with subsequent nu- 
clear capture ^^Fe {fi~ n) ^^Mn dominate. The time correlation of these neutrons largely 
reflects the capture rate of muons in iron (r = 206 ns) P8| . 

(iii) In the time interval < 6t < 20 ns there is an additional component, caused by 
muons which hit the outer veto and stop within the central detector. In this case, the time 
distribution corresponds to the muon time of flight from the veto to the central detector. 
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The solid histogram in figure ^ represents the expected time distribution from GEANT3.21 
simulations, which are in good agreement with the experimental data and are described in 
detail in E9| . 



Having identified events induced by cosmic ray interactions on iron using the outer veto, 
this background is strongly suppressed. The measurement indicated by full circles in figure |^ 
shows the remaining cosmic induced background, if sequences are rejected where the prompt 
events have simultaneously addressed modules in the central detector and in the outer veto 
system. These remaining sequences constitute the cosmic ray induced background for the Vf. 
analysis. They arise from the fraction of neutrons, which are produced outside the outer veto 
system, and are not absorbed in iron on their path to the detector (attenuation length of 
highly energetic neutrons in iron A = 21.6 cm[Q). The remaining spectrum consists of two 
components. The soft component is caused by neutrons from muon capture reactions and 
can be described as an exponential distribution e~^/^o with EqK,\A MeV. The much harder 
component is attributed to neutrons which have been produced in deep inelastic scattering 
processes of cosmic ray muons. This second component with a parameter of ^ 42 MeV 
covers the entire region of interest for the oscillation search. 

Compared to the background rates before the installation of the outer veto system (cor- 
responding to the energy spectrum with open circles in fig. |^), a background suppression by 
a factor 35 is achieved, resulting in a total rate of Res = (0.20 ± 0.01) mHz for the data 
cuts of the z/^— > Ve analysis in section |VT[ With this rate the cosmic induced background is 
smaller than the neutrino induced background. 

B. Neutrino induced background 

A second source of background reactions arises from the charged current (CC) and neutral 
current (NC) interactions of v^. and z/^ with the carbon nuclei of the liquid scintillator and 
iron nuclei of the inner passive shielding. To estimate the background contributions arising 
from different CC and NC reaction channels, the experiment takes advantage of having 
measured all relevant cross sections in a series of precision measurements [|l^, |19| . Thus, the 



calculated number of background events from conventional neutrino interactions does not 
rely on theoretical estimates of neutrino induced cross sections. This is especially important, 
as the i/-induced background is the dominant background contribution to the KARMEN 
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FIG. 9: Measurement of ""^^C ( z^e i e~ ) ^'^Ng.g, reactions (measuring points), leading to long lived 
coincidences between prompt e~ and delayed e"*" from ""^^Ng.s. decay (solid line (MC), shaded area 
(background)), (a) event time of e~, (b) visible energy of e~, (c) time difference between e~ and 
e+, (d) visible energy of e+. 

neutrino oscillation search. 

In the following we discuss the different //-induced background reactions in detail. For 
each background component we specify the experimental cross section as well as the detailed 
spectral information on energy and time, which have been used to calculate its contribution 
to the z/^— i> Ve oscillation search. 

I. The Vf, induced charged current reaction 

Exclusive charged current interactions of Vf. with ^^C can be detected by a delayed co- 
incidence consisting of a prompt electron from the inverse beta reaction ^^C ( z/g ,e~ ) ^^Ng.s. 
and the subsequent detection of a delayed positron from ^^Ng g, decays: 
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The lifetime of ^^Ng s. is r = 15.9 ms and the /3-decay endpoint amounts to Eq = 16.3 MeV. 
In total, 860 sequences of this type have been identified with a signal to background ratio 
of 61:1. Figure ^ shows the spectral information of the measured sequences, which are both 
for the prompt events and the delayed events in good agreement with the expectation from 
simulation. This fact underlines the reliability of the use of these simulated spectra in the 
likelihood analysis later applied. The measurements of KARMEN 1 and KARMEN 2 show 
full compatibility. For definiteness we use in the following the published CC event sample 
of KARMEN 1 , which leads to a cross section of 



It is the small fraction of 1.7% of ^^Ng s. decaying within the first 300 /is and depositing visible 
energies of less than 8 MeV which contribute to the expected background in the search. 
This background is extrapolated from the measured number of charged current sequences 
with time differences of 0.5 < At < 35.5 ms to the smaller time interval 0.5 < At < 300 /is on 
basis of the known ^^Ng.s, lifetime and the ^^Ng.s. energy spectrum. The uncertainties in the 
extrapolation correspond to 5% accuracy in the prediction of this background component. 

Charged current reactions of Uf. on iron with subsequent neutron evaporation from the 
excited iron nucleus ^^Fe ( z/g , e~ n) ^^Co have been investigated and simulated as possible 
background channel. Despite the rather high cross section calculation of a = 34.8 x 10~^^cm^ 
for this reaction channel and the significant number of target nuclei of the inner passive 
shielding Nt = 2.4 x 10^°, reactions on iron do not give rise to background in the t'^^ Ue 
analysis. The suppression of this channel is caused by the low efficiency of the electrons, 
which are produced with energies up to 35 MeV inside the steel, to be detected in coincidence 
with the neutron events inside the central detector. 

2. Random coincidences 

Neutrino-nucleus interactions as well as neutrino electron scattering increases the number 
of events in the positron time window in the first few /is after beam on target. This implies 
an enhanced rate of random coincidences between a neutrino induced (prompt) event and 



a = [9.4 ± 0.4(stat.) ± 0.8(sys.)] x lO'^^cm^ 



(11) 



23 



> 250 

CD 

- 200 

C/D 

I 150 
> 

CD 



50 




- J 


\ a) 


'- A 








- ^^^^^^^^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



400 



5 10 15 20 25 30 35 40 45 50 
energy prompt event [MeV] 




2 4 6 8 
visible energy [IVIeV] 



FIG. 10: (a) Measurement of neutrino induced reactions in the time window 0.6 < t < 
10.6 fis. The calculated contributions are (bottom to top): u — e'^ scattering, ^^Fe ( z^e ; e" ) ^^Co, 
( z^e , e" ) i^N*, {ue,e-) i^Ng.s.and ^^C ( i/ , jy' ) ^^c*. (b) Energy distribution of uncorrelated 
delayed events. 



a low energy event from environmental activity. Random coincidences, initiated by cosmic 
induced events, are accounted for in the measurement of the cosmic induced background. 
The probability P^c of an uncorrelated event to follow a neutrino-nucleus interaction as well 
as its spectral information are extracted by applying the search criteria to uncorrelated 
events, for example to events recorded in preceding beam periods. This method allows to 
determine Pre and the spectral information of the delayed events with high statistics. 

In the energy range from 11 < < 50 MeV and in the time window 0.6 < t < 10.6 fis 
after beam on target, 1567 neutrino-nucleus interactions are measured. The neutrino inter- 
actions arise mainly from two different types of neutrino-nucleus interactions. The largest 
contribution arises from the inclusive charged current reaction ^^C {Uf. ,e~ ) ^^N, as well as 
from neutral current reaction ^"^C {u ,1/' )^'^C* (1"^,1;15.1 MeV) with u = {vei'i^ii)- The 
neutral and charged current contribution are clearly visible in the energy spectrum of the 
measured neutrino- nucleus interactions (see Fig. |T0|). The delayed events of random coinci- 
dences are uniformly distributed in time, and their energies are close to the threshold of a 
single detector module (mean energy {E) = 1.1 MeV) as shown in figure |lO|(b). 

The probability P^c for an uncorrelated event to occur with a time difference of up to 
5 < At < 300 /is and within a coincidence volume of = 1.3m^ after the prompt event is 
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determined to be Pre = (5.5 ± 0.4) x 10 ^ |]65[. 



The expectation value for the neutrino-induced random background Nrc is obtained by 
multiplying the number of measured neutrino-induced reactions N^, with the probability Pre- 
Using this method, the statistical accuracy of Nrc is 7%. The measured spectral information 
is used for the likelihood analysis. 

3. i>e contamination 

The only background source, which can not be directly extracted from the data, is the 
contamination of the neutrino beam with Ve produced in the vr^- yU~ decay chain. Detailed 
Monte Carlo simulations, including a three-dimensional model of the ISIS target, and its 
surroundings are used to obtain the fraction of tt^ and /i^ decaying before they undergo 
capture on nuclei of the target materials [^, |52|. The overall ratio of produced in the 



ISIS target relative to from /i^-decay amounts io e = 6.4 x 10"'^. This ratio is further 
reduced by taking into account, that the lifetime of /i~ depends on the target material and 
is in general shorter than the /i^ decay time (see figure [^), leading to a further reduction of 
z/g by a factor of 0.764 in the time window of 0.6 < t < 10.6 yus. Finally, the Ve spectrum 
from /i~ decay (Fig. |l](b)) leads to a lower flux averaged cross section of a = 72.0 x 10~^^cm^ 
for the p ( z/g , e"*" ) n reaction (see table |). Taking all effects into account, the intrinsic z/g 
contamination leads to the smallest background contribution in the z/g search. 

C. Beam correlated neutron background 

Each 800 MeV protons of the ISIS beam produces typically 25 spallation neutrons in the 
target with energies up to 400 MeV The 7 m steel shielding between ISIS target and 
detector reduces the neutron flux by a factor of more than 10^^. Despite the flux reduction, 
punch-through neutrons are observed in the central detector. However, these high energy 



neutrons closely follow the ISIS double proton pulses |20[ and are restricted to the time 
window of t < 500 ns after beam on target. Setting the lower time cut for the positron 
window at tpr, > 600 ns after beam on target, completely eliminates reactions from these 
neutrons. 
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VI. DATA REDUCTION 



A. Raw data 

The results presented here are based on measurements from February 1997 to March 2001. 
During this time, protons equivalent to an accumulated total charge of 9425 Coulombs have 
been stopped in the ISIS target. This corresponds to a total number of 

= 2.71 X 10^^ (12) 

neutrinos for each of the flavors z/g, i/^ and z/^ produced at the ISIS beam stop. 

In total, the KARMEN data acquisition system recorded 3.7 x 10^ events. Out of these 
single events, 1.93 x 10'' have no hits in the veto counter system and deposit more than 
11 MeV and hence can be classified as candidates for a prompt event of a delayed coincidence. 
Requiring in addition the detection of a second event without veto hits in the following 500 /is 
results in 3.5 x 10^ delayed coincidences. After application of the general event requirements, 
defined in section the sample size shrinks to 3464 coincidences with more than 99% of 
these coincidences outside the time window of the z/^— analysis. 

The detector system was 777.4 days online, excluding additional measurements for specific 
background studies and calibration purposes. Taking into account ISIS beam on times, the 
duty cycle, and a 10 /is long neutrino time window, the effective neutrino measuring time 
amounts to 7.5 hours. 



B. Final selection criteria 

The final selection criteria have been evaluated in order to optimize the sensitivity of the 
experiment. Since the true values of the oscillation parameters are unknown, we optimized 
the data reduction to deliver the most stringent upper limit on sin^(26) for a given Am^ 
under the assumption that there are no t'^— Ue oscillations. Even a small oscillation signal 
would then first materialize as a much less stringent upper limit then the experimental 
sensitivity. The optimized cuts were obtained by simulating and analyzing experimental 
outcomes with different cuts leading to different event statistics It turned out that the 



achievable sensitivity only slightly depends on the variation of reasonable data cuts. 
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The final data cuts are as follows: Accounting for the ISIS time structure, the e^ from 
i^fi—^ oscillations must be detected in the time interval of 0.6 < tpr, < 10.6 /is after beam on 
target, in which 84.0% of all z/^ are expected. The lower time cut of 600 nanoseconds is chosen 
to eliminate any contributions from beam correlated fast neutrons (see section |V C|) . The 



lower cut on the visible energy deposit Ep^. of a positron candidate is 16 MeV. This energy 
cut eliminates the neutral current contributions ^^C f z/ , z/' ) ^^C* to the neutrino induced 



random background (fig. |T0D and also suppresses the soft component of the cosmic induced 
background (fig. 0). No fiducial volume cut for the e"*" is applied. 

The time difference for the delayed neutron capture event is restricted to the interval 
5 < At < 300 /is. Here, the lower time cut is fixed by a minimum hardware deadtime 
after the electronic read-out of the prompt event. The upper time cut at At < 300 fis 
is an outcome of the MC procedure mentioned previously and reflects the different time 
distributions of delayed events from neutron capture (r ^ 120 fis) and from the background 
reactions of random coincidences (uniformly distributed) and charged current coincidences 
(r = 15.9 ms). 

The remaining data cuts for neutron capture events are the coincidence volume of Vc = 
1.3 m^ and a maximum energy of the neutron capture event of Edei. < 8.0 MeV. Table |I| 
gives a summary of the applied data cuts and the corresponding efficiencies e, resulting in 
a total efficiency 

Stotii^,) = 0.192 ± 0.0145 (13) 
for an oscillation signal at large Am^. 

C. Data reduction 

Applying the final selection criteria to the entire KARMEN 2 data set results in 15 
candidate events. 

The total background expectation amounts to N'j^q = (15.8 ± 0.5) events for the compo- 
nents described in section 0. As can be seen from the summarizing table III , the background 



is dominated by neutrino induced processes, whereas the cosmic induced background con- 
tributes to only 25% of the total rate. The relative uncertainty of the background expectation 
amounts to 5%, reflecting the accuracy of the in situ measurement of the three dominating 
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event 



data cut 



efficiency e 



check on 
previous history, 



(see sec. [TVl) 

0.6 /US < tpr. < 10.6 /US 



16 MeV < Epr. < 50 MeV 



0.709 

0.840 
0.775 



("-,7) 



5 /US < At < 300 /is 
Edei. < 8.0 MeV 
Vr = 1.3 



0.416 



TABLE II: Final data cuts and efficiencies for the f^— > fe search. The efficiency for the energy cut 
corresponds to oscillation parameters Am^ > 100 eV^. 



background 


expectation Ni 


method of determination 


Cosmic induced background 


3.9 ±0.2 


measured in diff. time window 


Charged current coincidences 


5.1 ±0.2 


measured in diff. energy, time windows 


Ve ind. random coincidences 


4.8 ±0.3 


measured in diff. time window 


i/f. contamination 


2.0 ±0.2 


MC- simulation 


Total background N'j^q 


15.8 ±0.5 





TABLE III: Expected background contributions 



background components in different energy and time windows. Figure |Tl| shows the spec- 
tral distribution of the 15 candidate events with the superimposed background expectation, 
normahzed to 15.8 events. In each plot the measured data agree well with the expected back- 
ground distributions. There are no obvious deviations from the background expectations, 
neither for the prompt nor delayed events. 

Already, the agreement of the number for measured events with the expected background 
does not give any hint for an oscillation signal within the KARMEN 2 data. In the following, 
we will set upper limits on the oscillation parameters, also using spectral information of the 
candidate events. 
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FIG. 11: Final event ensemble (a) time of prompt events, (b) energy of prompt events, (c) time 
difference between prompt and delayed event, (d) energy of delayed events, (e) spatial correlation 
and (f) distance to target of prompt event. The 15 oscillation candidates are in very good agreement 
with the background expectation of 15.8 events (solid line). 



VII. DATA ANALYSIS 

For z/^^ z/e oscillations with maximal mixing (sin^(20) = l) and large mass differences 
(Am^ > lOOeV^), an oscillation signal of (2913 ±269) sequences is expected (see table [TVl ). 
This number includes a small contribution from z/^ produced at the intermediate ISIS /xSR 
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detection reaction 



expectation N^ig neutrino source 



p ( , e+ ) n 
p ( , e+ ) n 
i2C(Pe,e+ n) 



2716 ± 268 
73 ±7 
125 ± 17 



P^— > from main target 
i^/x— > i^e from ^SR target 
i^fi^ from main target 



Total iV,7/(sin2(2e) = 1, Am^ = lOOeV^) 



2913 ± 269 



TABLE IV: Expected f^— > Vf. oscillation signal for maximal mixing 



target [^. The systematic error of the oscillation expectation is dominated by the neutrino 



flux uncertainty of 6.5% |52| and the error in the determination of the neutron detection 
efficiency of 7.0 %. 

Having measured 15 events with a background expectation of 15.8 events, there is no 
indication for the presence of an oscillation signal in the KARMEN data. Ignoring, in a first 
step, the spectroscopic information of the measurement and interpreting the experimental 
outcome as a pure counting experiment, an oscillation signal larger than A^^ig = 7.4 events 
is excluded in 90% confidence interval (C.I.) [0, However, such a simplified approach 
does not make any use of the spectroscopic quality of the data. In order to extract more 
information on a potentially small oscillation signal in the final event ensemble, a single 
event based maximum likelihood analysis is applied to this ensemble. 



A. Likelihood analysis 

The purpose of a maximum likelihood analysis is the separation of a potential signal from 
background by maximizing the likelihood function with regard to some unknown parameters. 
In this case, the signal corresponds to (e"'",n) sequences from z/^— Ve oscillations, the unknown 
estimators are the oscillation parameters sin^(20) and l\rv? . The likelihood function L is 
defined as: 

^ sam'ple. 

L{r,^m^)= J] [r-/,,,(f„,Am2) + (l-r)-/bg(fO] (14) 

71=1 

with the following definitions: 

• The event sample with A^samp/e=15 candidate events is characterized by the information 
on the energy E^r and time Tp^ of the prompt event, the energy of the delayed event 
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Edei and the time difference AT and position correlation Ax between prompt and 
delayed event. This information for each event sequence n is represented by the vector 

•^n (-^pr; Tpj., Efi^i, AT, Ax). 

fsig and fbg are the probability density functions for the vector Xn in case of event n 
being a signal or a background event. 

The parameter r describes the signal fraction in the data and is connected to sin^(20) 
by the linear transformation 

«in^(2e) = ^e.p., A 2, (15) 

^si/(sin (26) = l,Am') 

with the calculated oscillation signal N^^^' {sin^ {2Q)=l,Am'^) for maximal mixing as 

shown in figure (a). 

Assuming no correlation for the j = 5 observables of x, the probability density function 
is factorized to: 

5 

i=i 

= f{E,r, Am2) . f{Tpr) ■ f{Edel) " /(AT) ■ /(Ax) 

Due to the small event sample size of 15 events, the fit is not performed by varying 
simultaneously the signal and all background components individually. In contrast, 
the four individual background components are added up to one total background 
component 

/.,=x^Q- (n/,,,j (17) 

with the coefficients q being the expected relative contributions of the background 
channels. The values of q are given by the ratio of the expected number of background 
events Ni of each component and the total background expectation N^^q (see table 
pni ): Ci = Ni/N'^Q, thereby satisfying the normalization condition J^i^^i = 1- 

With the normalization constraint of the probability density function, the parameter 
r determines also the background contribution in the likelihood maximization: 

Nbg = (1 - r) ■ Nsample (18) 



31 



with the above definitions, maximizing L with regard to sin^(20) and Am^ is a pure shape 
analysis and does not take into account the knowledge of the total background expectation 
iV^o'. To include this quantitative information, the likelihood function is weighted with 
a Poisson probability term Pp computing the probability of measuring Nhg{r) background 

rexp 

''bg 



events for an expectation of N'j^J! events: 



Cir, Am2) = £(r, Am^) • Pp (iV,,(r)|iV,7^) (19) 

with 

The expansion in the Poisson probability from the discrete factorial n\ to the Gamma func- 
tion r(a;) = Jq°° e~^t^~^dt with T{n + 1) = n\ ensures a continuous calculation of the Poisson 
probability for any signal ratio r. 

Maximizing the above defined likelihood function C for the final KARMEN 2 data results 
in a best fit for r compatible with the no-oscillation solution. In fact, the global maximum 
of C is reached slightly in the unphysical region, at oscillation parameters: 

sin2(2e) = -2.4 ■ 10~^ Am^ = 5.4 eV^ (21) 

Restricting the analysis to the allowed region, the likelihood function is found to be maximal 
at: 

sin2(2e) = 8.0 • 10-^ Am^ = T.OeV^ (22) 

Table |V] shows the number of (e"'",n) sequences from z/^— > z/g oscillations of some selected 
parameter combinations k. Also given are the differences of the likelihood values to the 
maximum in the physically allowed region: — Aln£ = InCk — In Cmax- The logarithmic 
likelihood value of the best fit differs from the likelihood value for no oscillation by only 



0.21 units. As will be discussed in section VII B, a statistical analysis of the likelihood 



function indicates that for boundaries of 90% confidence intervals (C.I.), typical differences 
of — Aln£ 4 — 5 have to be applied. This underlines the fact that the maximum at 
(sin^(20), Am^)max = (8.0 ■ 10""^, 7.0 eV^) is statistically in excellent agreement with the 
null hypothesis of no oscillations. 

Furthermore, simulations of comparable event ensembles, with no oscillation signal but 
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sm2(2G) 



-AlnC 



no oscillation 0.21 

0.1 -1.7-10-2 -0.3 0.21 

5.4 -2.4-10-3 -4.4 -0.75 

7.0 8.0 - 10-^ 1.5 0.00 

100 2.1 - 10-^ -0.6 0.18 



TABLE V: Signal event numbers for selected oscillation scenarios. The values of — Aln£ indicate 
the difference of the likelihood value to the maximum in the physically allowed region (see text). 

background only, show that a global maximum at slightly unphysical oscillation parameters 
as it is in case here (Eq.|21]) is a typical result of the likelihood analysis of small event samples 

{-^sample) 15.8. 

Figure 0(a) shows the expected oscillation event numbers as a function of Am^ for 



maximal mixing sm^{2Q) = 1. In contrast, figure |T2|(b) demonstrates the results of the 
maximum likelihood analysis. For 90 slices per decade in Am^, the number of oscillation 
events for maximal likelihood N^^^ is plotted (solid histogram). For low as well as high 
values of Am^, the corresponding best fits are almost identical with the physical boundary, 
with values of N^^^ = —0.3 and N^^^ = —0.6, respectively. In a region of about 3 < 
Am^ < 30 eV^, stronger variations of the energy spectrum of a potential signal come into 
play: Since KARMEN has an excellent energy resolution of ~ 2 % for positrons with E ^ 
30 MeV, statistical variations in Epr of the small event sample can be easily interpreted by 
the likelihood analysis as modification of the background energy spectrum due to oscillations 
with an oscillation length comparable to the distance target-detector: 

Lose = 1 — o ~ 17m (26) 

For energies 12 < E{u^) < 52.8 MeV, equation ( P3[ ) leads to oscillation parameters of about 



3 < Am^ < 15eV2||U^. It is important to note that the results given in figure |T^(b) for 
3 < Am^ < 30 eV^ are statistically perfectly compatible with the no-oscillation solution, as 
will be discussed in the next section. 
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FIG. 12: (a) Expected oscillation signal for maximal mixing sin^(20) = 1. (b) Results of the 
likelihood analysis: The solid line shows the best fit of a Vp.^Vf. signal in the KARMEN2 data. 
The dashed line corresponds to the upper bound of the derived 90% confidence interval (C.I.) for 
an oscillation signal. Note that there is no lower bound of the 90% C.I. for all Am^. 



B. Upper limits on oscillation parameters 

Finally, the confidence intervals for the parameters sin^(26) and Am^ have to be deduced 
from the experimental likelihood function. Recently, there have been discussions about 



various approaches in order to obtain confidence regions, especially under the aspects of event 
samples of low statistics, oscillatory behavior of the likelihood function as well as parameter 
determination near physical boundaries. In the following, we adopt the Unified Approach 
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|57| which is a frequentist approach with a specific ordering principle: In the [sin^(20), Am^]- 
plane, a 2- dimensional confidence interval (C.I) for the oscillation parameters is constructed 
from the comparison of the experimental likelihood value Aln£ = In £(sin^(20), Am^) — 
In £((sin^(20), Am^)ma^. with the outcome of a large sample of Monte Carlo simulations 
of so-called toy experiments for this term. These simulations are based on the detailed 
knowledge of all resolution functions and the spectral information on the background. In 
addition, they comprise the expected experimental signal for an oscillation hypothesis with 
given parameters (sin^(2B), Am^). The hypothesis is then accepted in the 90% C.I. if 
the experimental value does not lie within the outer 10% tail of the simulated —A In re- 
distribution. For a complete statistical analysis, the entire parameter space [sin^(20), Am^] 
is scanned to extract the according region of confidence. 

In figure |l2|(b), the result of this approach is shown in terms of excluded oscillation events. 
The dashed line corresponds to the limit of the 90 % confidence interval, excluding larger 
signal event numbers. For Am^ = lOOeV^ an oscillation signal stronger than 5.1 events is 
excluded in the 90% C.I., while for low Am^< 0.1 an oscillation signal larger than 6.0 events 
is excluded. Though one of the major features of the Unified Approach is the possibility 
of extracting lower limits within the same analysis, no lower limit of the 90 % C.I. appears, 
demonstrating the compatibility of the likelihood result with the no-oscillation hypothesis 
regardless of the chosen value for Am^. 

The exclusion plot in the 2-dimensional [sin^(2B), Am^]-plane (fig. |T^)is derived by di- 
viding, for all values of Am^, the excluded events (see the solid line in figure [l^ (b)) by the 
expectation for maximal mixing (figure [T^ (a)). This results in the 90% C.I. limits: 

sin2(20) < 1.7- 10"^ Am2>100eV2 (24) 
Am2 < 0.055 eV2 sin2(2e) = 1 (25) 



The entire exclusion curve is shown in figure |I3| as solid line, excluding parameter combina- 
tions in the area right to the curve. 

An important criterion of an experimental result and a derived upper limit is the question 
of how close the limit quoted is to the experimental sensitivity. Following [^], the sensitivity 
is defined as expectation value for the upper limit of the 90 % confidence interval under the 
assumption of no oscillations and is gained by simulations of experiments' outcomes. The 



KARMEN2 sensitivity as a function of Am is shown in figure |I3| as dashed line. The 
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FIG. 13: KARMEN2 90% C.I. result deduced with the Unified Approach (solid); 90% C.I. sen- 
sitivity within the Unified Approach (dashed) and 90% C.I. in the Bayesian approach (dotted). 
Regions to the right of the curves are excluded. Note the zoom of the axis in sin^(20), not reaching 
up to 1. 

sensitivity (sin^(20)) for Am^ = 100 eV^ amounts to 

(sin2(2e)) = 1.6 • IQ-^ 90 % C.I. (26) 

For completeness, we also perform a Bayesian approach to derive an upper limit on the 
oscillation parameters sin^(20) and Am^. In the Bayesian framework, the upper limits for 
fixed Am^ are obtained by integrating the likelihood function L. This integration implies 
the use of a prior probability density distribution for sin^(29) and decomposes the 
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2- dimensional problem into a one-dimensional treatment. We used a uniform prior in a log- 
arithmic metric of the oscillation parameter sin^(2B). In both the frequentist and Bayesian 
approaches, we restrict the parameter space to the physically allowed region. The Bayesian 
90% C.l. approach yields more stringent upper limits shown as dotted line in figure [T3| with 

sin2(2e) < 1.3 ■ 10"^ Am^ > 100 eV^ (27) 

Because of the ambiguities in choosing the probability density distribution for sin^(2B) 
as well as the 2-dimensional oscillatory behavior of the likelihood function, we do not favor 
the Bayesian extraction of confidence intervals but refer to the results deduced within the 
frequentist Unified Approach (see Eq. ^). The resemblance of the KARMEN exclusion 
curve with its sensitivity curve underlines the fact, that the likelihood analysis results in no 
indication of a z/^— > Ve oscillation signal in the KARMEN 2 data. 



C. Comparison with LSND and other experiments 

The parameter space for oscillations z/„— > z/e excluded at 90% C.l. by the KARMEN 2 



measurement is shown in figure |TJ. The KARMEN result sets the most sensitive limits so 
far on t'^— > Ve oscillations in the parameter region of 0.3 < Am^ < 30 eV^. At higher Am^ 
values, the area right to the right to the exclusion curve is also excluded by a combined 



^ii^'^e and Vfi,-^Ve search of CCFR The search for z/g disappearance at the Bugey 

reactor ||60| excludes small Am^ values, at large amplitudes A > 0.03. The parameter 
area excluded by KARMEN covers large parts of the parameter combinations favored by 
the LSND experiment The LSND result plotted here shows areas obtained by cutting 
the experiment's logarithmic likelihood function at constant values 2.3 and 4.6 units below 
the likelihood maximum |]I2[. For values of Am^ < 2eV^, the oscillation signal expected 



in KARMEN based on the LSND region (In C^ax — 2.3) corresponds to a range of 3 to 
14 oscillation events. As shown in figure |12|, a signal larger than 6 events is excluded at 
90% CL . At Am^ > 20 eV^, the expected LSND signal of 7 to 13 oscillation events in 
KARMEN is in clear contradiction to the KARMEN upper limit of 5.1 (6.5) events at 90% 
C.L (95% C.L). 

These examples based on expected additional Ue events from z/^— z/g demonstrate that at 
smaller values of Am^ there is a restricted parameter region statistically compatible with 
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FIG. 14: Comparison of oscillation searches performed by different short basehne experiments. 

both experimental results. At high Am^ values, the LSND solutions are in clear contradiction 
with the KARMEN upper limit. 

VIII. CONCLUSION 

Results based on the entire KARMEN2 data set collected from 1997 through to 2001 
have been presented. The extracted candidate events for Vf, are in excellent agreement with 
background expectations showing no signal for v^-^ z/g oscillations. A detailed likelihood 
analysis of the data leads to upper limits on the oscillation parameters sin^(2G) and Am^ 
excluding parameter regions not explored analyzed by other experiments. 
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These limits exclude large regions of the parameter area favored by the LSND experi- 
ment. A more quantitative statistical statement on the compatibility between KARMEN 
and LSND has to be based on a combined statistical analysis of both likelihood functions 
pT|. Such a detailed joint statistical analysis has been performed 



The negative search for z/g from muon decay at rest presented here sets also stringent limits 
on other potential processes of Ve production such as lepton family number violating decays 
yU"'' — > e"^ + z/e + or neutrino oscillations z/g— > which will be discussed in a separate 
paper. Future experiments such as the MiniBooNE experiment at Fermilab aim at 
investigating the LSND evidence and the oscillation parameters not yet excluded by the 
z/^— z/e search presented here. 
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